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ABSTRACT 

A large-scale closed magnetic field can transfer angular momentum and energy 
between a black hole (BH) and its surrounding accretion flow. We investigate the 
effects of this magnetic coupling (MC) process on the dynamics of a hot accretion 
flow (e.g., an advection dominated accretion flow, hereafter ADAF). The energy 
and angular momentum fluxes transported by the magnetic field are derived by an 
equivalent circuit approach. For a rapidly rotating BH, it is found that the radial 
velocity and the electron temperature of the accretion flow decrease, whereas the 
ion temperature and the surface density increase. The significance of the MC 
effects depends on the value of the viscous parameter a. The effects are obvious 
for a = 0.3 but nearly ignorable for a = 0.1. For a BH with specific angular 
momentum, a^, = 0.9, and a = 0.3, we find that for reasonable parameters the 
radiative efficiency of a hot accretion flow can be increased by ~ 30%. 
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Introduction 



As a variant of Blandford-Znajek (BZ) process (IBlandford &: Znajekl 119771 ). the mag- 



tion disks has received much attention (e.g. 



Blandford 


1999; 


Li & Paczvnski 


2000: 


Li 


2002 
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Wang et al. I l2002h . By virtue of the large-scale closed magnetic field lines that connect the BH 



the BH and the disk. ( 


2002) 


1973; 


Novikov & Thorne 


1973; 



Page fc Thorndll974j ) case the MC process may change the 



local radiative flux significantly. 



Apart from SSD, ADAF is anot her important model of the accretion flow (jNarayan fc Yi 



1998 



1994, 1995 ; Abramowicz et al1ll995 ; see Narayan. Mahadevan &: QuataertlEgga and lKato. Fukue fc Minesh 



for reviews). It has been apphed to a number of accreting BH systems and successfully 
explains their spectral characteristics (see Narayan 2005 and Yuan 2007 for recent reviews). 
In ADAF models the thickness of the accretion flow is of the same order as ra dius, i.e. H ^ r, 
which means the large- scale field is easier to form in an ADAF than in a SSD (ITout fc Pringk 
19961 ; iLivio et al.lll999l ). So it is interesting to investigate the influences of the MC process 
on an ADAF. Very recently, lYe et al.l (120071 ) discussed this problem based on the self-similar 
solution of the ADAF. In this paper, we investigate the influences of the MC process on 
ADAFs through global solutions. 

To properly assess the dynamical effects of large-scale magnetic fields on the accre- 
tion flow, one needs to obtain the fluids and the fields at the same time by solving the 



transfield equation, which is a non trivial non . 

unct i ons JUzdenskv 



gular surfaces a nd free 
simulations 



Haw 



McKinnev fc Gammie 



e^ 



200 



2000 

T 



inear parti al differential equation with sin- 



Hawlev &: Bal 



Hirose et al. 



DU 



2004 

T 



2002 



20051). An alternative way is MHD 



Koide 2003; De Vilhers et al. 2003 



20041 ). However, both of these approaches are com- 



plicated. Lai ( 19981 ) and Le 3 (Il999al . Ibh adopted a phenomenological approach to research 
the magnetic coupling between the neutron star and its surrounding slim disk. They speci- 
fied an ansatz for the magnetic fields and then numerically solve the basic equations of the 
accretion flow. In their model the disk is geometrically thin, 8/ Or ~ 1/r ^ 1/H ~ d/dz, 
and so the expressions of electromagnetic forces can be much reduced by omitting d/dr 
terms. However, an MCADAF is thick and the expression of the electromagnetic force is 
complex. Here for simplicity we treat the MC process as a source of energy and angular 
moment without considering the radial and vertical components of the electromagnetic force 
in the momentum equations. 

We derive the energy and angular momentum fluxes i n the Kerr metric by using the 
approach of equivalent circuit (IMacdonald &: Thornd 119821). But for si r nplicit y a pseudo- 
Newtonian potential of a rotating black hole given by iMukhopadhyayl (120021 ) is adopted 
when we solve the solutions of the accretion flow. 



In Section 2, we describe the MCADAF model and calculate the energy and angular 
momentum fluxes transferred by the magnetic fleld. In Section 3, we write down the basic 



3 



equations describing the MCADAF. The numerical results are presented in Section 4 and 
Section 5 is devoted to a summary and discussions. Throughout this paper the geometric 
units c=G=l are used. 



2. MCADAF Model 

We assume the ADAF is stationary and axisymmetric. The ADAF extends from the 
outer edge, Vout, to the BH horizon rn- There are two kinds of magnetic fields in this model, 
i.e., large-scale closed magnetic field that connects the BH with the ADAF and small-scale 
tangled magnetic field, with the former contributing to the MC process and the latter to 
the viscosity. We assume these two kinds of fields work independently. If not mentioned we 
refer the magnetic field as the large-scale closed one hereafter. The region between the BH 
and the ADAF is assumed to be ideally conducting and force-free. 

The field lines are supposed to distribute in the ranges of {rH,rout) on the disk and 
(0, ^o) on the horizon. Due to the lack of knowledge about the magnetic field around the 
BH, we assume that the field threading the BH is constant, i.e. Bh{0) = const. The field 
threading the ADAF is assumed to decrease with r following a power law form, but within 
the marginally stable orbit, < r^s, the radial velocity of the accretion flow increases much 
faster thus the filed is likely to increase with radius. Given this consideration, we assume 
the field has the following distribution, 

n (r■^- R F(r\ - / ^0 exp(r/rp - 1) for r^^ < r < 

Here rH = M{l + y/f — al^ denotes the radius of the BH horizon, a* is the dimensionless 
spin parameter of the BH, and rp = th + A(rms — th)- 



Moderski. Sikora fc Lasotal (119971 ) gave an estimation of Bh with the balance between 
the ram pressure of the falling material and the magnetic pressure, i.e., Bfj/Sn r-u p ^ 
Mr>/ (47rr|^). Since the ram pressure can be larger than the magnetic press, we introduce a 
parameter cb to indicate the strength of the magnetic field threading the horizon as 

Bh = CBV2Kl/rH, < cb < 1. (2) 



In the following derivation of this subsection, the Boyer-Lindquist coordinates are used. 
Assume all the field lines threading the BH are connected with the disk, then from the 
conservation of magnetic flux we have 
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^ = y" BH{pvo)Hded<i) = Jb, (^^^ drd^, (3) 

where the subscripts "H" and "D" are used to indicate the quantities on the horizon and 
the equatorial plane of the disk {6 = vr/2), respectively. The Boyer-Lindquist coordinates 
are given as 

j:^ = {r^ + alMy -alM^ A sin^e, = + ajM^ cos^ 6, 

A = + alM^ - 2Mr, w=(S/p)sin^. ^' 

Since A = at r = r/^, the lower boundary of the integration interval in the second equality 
is set to be rn + Sr, where 6r is a small quantity and taken as 5r = 0.010. Substituting 
equation ([T]) into equation (JSj) we get 

J BH{pzu)Hded(j) 2MrH{l-cosdo)BH ^ ... 
Bo = J r = J r = 2rHk{a^,n)BH/M. (5) 

/ Fir) (^) ^ drd<l> J Fir) (^) ^ rfrrf0 
Given the configuration of the field, we can derive the energy and angular moraentum 



flux in the MC process by using the modified equivalent circuit approach (jWang et al.ll2002l ). 
Considering a loop corresponds to two adjacent flux surfaces (characterized by the magnetic 
flux and + A\&), the electromotive force due to the rotation of the BH and the disk are 
expressed as 



Aeh = iA'$/2'K)QH, Ae 



D 



iA'^/2n)Q, A^i = 2Tiiwp)HAe ■ Bh- 



(6) 



The minus sign in the expression of Ae£, arise from the direction of the flux. The parameter 
VL is the angular velocity of the ADAF, VL^ = a^:/i2rH) is the angular velocity of the BH 
horizon. 



The equivaleii t surface resistivity of the BH horizon is An ( IMacdonald fc Thorndll982 
Thorne et al.lll986l ). while the surface resistivity of the disk is ~ 1/ iHa) = 4:7rr)/H, where r) = 



l/(4:7ia) is the diffusivity of the magn etic field. As in many papers (e.g. iLubow et al.lll994j : 
Lovelace et al.lll995l : ISoria et al.lll997l ). we assume rj to be of the same order as the Shakura- 
Sunyaev (1973) kinematic a-viscosity coefficient, i.e., r/ ~ z/ = aCgH. The resistances of the 



^The influence of Sr can be ascribed to cb as tfie effects of tlie MC process are mainly determined by the 
strength of the field in the region r > Vp (we will show this in Sec. 4). 
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annular ring on the horizon and the disk are thus 



AZ„^4..^^^H^^. (7) 



1 Ar 2aCs ■ Ar 

Az£) 



Ha 27izu£) Wo 
Thus the current in the loop is 



J _ AeH+Aep _ /'A*\ Qr-^ _ _!_ a*{l~PHD) j\ fr> 

AZh+AZd y 27T ) AZh-{1+0 1+€ ■ 2csc2 0-l+v^]^ ' 
(<: = ^-^g = acsTAjjj I Ar I . <q _ _Q_\ 

— AZh Ph'^d I Ae I ' t^HD — Qh ' 



(9) 



In order to obtain ^ and / we have to find the value of Ar/A6' , which is related to the 
mapping relation between the angular coordinate on the horizon and the radial coordinate 
on the disk, i.e. 9{r). According to the conservation of the magnetic fiux between the 
adjacent two fiux surfaces, 

dm = Bh- 2-K{wp)^de = -B, ■ 27T (wp/VA^^dr. (10) 
Substitute equations ([H) and (I5l) into the above equation, we get 



M — ; — = kia^.n)- , Fir) = G[a^,r,n). (11) 

Integrate equation (fTTj) and we obtain the mapping relation: 

cos ^ = cos ^0+ / G{a^:,r,n) ■ dr. (12) 

JrH 

Our calculations show that, for the ADAF, the ratio of the height to radius around is 
< 0.3, thus we assume 9o = OAtt so that cotOo ^ 0.3. Substitute equation (ITTl) into equation 
([9]) we have 

2aCsG{a^:,r, n)~^ 



2 csc2 ^ - 1 + ,/l - a? 



(13) 



Since the current / on the BH horizon feels Ampere's force, the BH exerts a net torque 
on the magnetic fiux tube 

^Tmc = tuBhIpAO = — — / = y . ■ MBfj ■ Ar. (14) 



27r 



1 + e) (2 csc2 6-1 + v^r^^ 



-6- 



Prom the second equality it is easy to find that this torque equals to the torque exerted 
on the disk by the same flux tube, or equivalently speaking, the angular momentum flows 
between the BH and the disk. The angular momentum flux can be written as 

The power transmitted to the disk through the tube is given by 



APmc = /AsD + J^AZij = 47rr//A^c^^Ar + 47rr//Mc(^^F-^^)Ar = I^PMw + ^Qohm, (16) 
where 

- AZn + AZn ' ^^^^ 
is the angular velocity of the magnetic flled hues, /S.Pmw = 'inrHMc^^r is the rate 
of the mechanical work done by the electromagnetic torque on the disk and AQohm = 
AnrHMci^F — f^)Ar is the rate of Ohmic heating in the disk. 

It is easy to calculate the power dissipated on the BH's stretched horizon intersecting 
with the flux tube: 

AQbh = I'AZh, (18) 
which increases the irreducible mass of the BH. 



3. Basic Equations of the Accretion Flow 

We assume the energy and angular momentum transferred by the MC process deposit 
into the accretion flow homogeneously in the vertical direction. The height- averaged basic 
equations describing the MCADAF can be written as 

M = -inrpHv = const, (19) 

M-^{nr^) + ^nrHMC = {^7rr\^H) , (21) 
pvT.^-{l-S){qL + ^)-<l^e, (22) 
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pvTe 



dSe 

dr 



Q 



Ohm - 

2ir' 



Qie - q 



(23) 



Here M is the accretion rate, Cg = \ppTp is the isothermal sound speed, p = Vgasj Pt = pc^/Pt 
is the total pressure of the tangled magnetic field and the gas pressure, Pt is the ratio of 
the gas pressure to the total pressure and is fixed at its "typical" value Pt = 0.9, T is the 
temperature, s is the entropy. The subscripts "i" and "e" indicate the quantities for ions and 
electrons, respectively. The quantit y Try, = —ap is the rjp com ponent of the viscous stress 
tensor adopting the a prescription (IShakura fc SunyaevI Il973l ). H = Cs/flx is the vertical 
scale height, Qk i s the Keplerian a ngula r velocity calculated by using the pseudo- Newton 
potential given by iMukhopadhyayl (l2002l ). 6 describes the fraction of the total energy that 
directly heats the electrons and is set to be 5 = 0.3 following the detailed modeling r esult to 
the supermassive black hole in our Galactic center (lYuan. Quataert fc Narayanll2003l ). = 
rTry, [dfl/dr) is the heating rate of the viscosity, qie represents the volume energy transfer 
rate from ions to electrons via Coulomb collisions, q~ is the cooling rate of the electrons , 



which consists of bre msstrahlung, synchrotron, and Comptonization (INarayan &: Yil Il995 
Manmoto et al.l 119971 ) . and Qohm = ^Qohm/4:7!'rAr = Hmc{^f — ^) is the rate of Ohmic 
dissipation per unit of the disk. 

Adopting the no-torque boundary condition at the horizon, we integrate equation 
from rn to r and get the conservation equation for the angular momentum 

(r)=No, (24) 



arc; 1 
I + + -^Tl^c\ 



where 



V 



TMcir) = 



dTi 



MC 



dr 



dr, 



Nn 



lo + -^Tljcirn) = const, 
M 



(25) 
(26) 



with Iq being the angular momentum per unit mass swallowed by the BH. The three terms 
on left-hand side of equation correspond to the advected angular momentum, viscous 
torque and magnetic torque due to the field lines in the range from r to rout- 



One more relation is given by the equation of state, 

Pgas = k(Ti/fii + Te//ie)/m^ 



(27) 



where p is the mean molecular weight, k is the Boltzmann's constant, and m„ is the atomic 
mass unit. 

Thus we have a set of six equations including one integral, two algebraic, and three 
differential equations, i.e. equations (|T9l) . fl20l) . fl22l) - fl2^ . and fl271) for six unknown quan- 
tities, H (or Cs), V, p, Q, Ti, and Tg. This set of equations can be solved with three outer 
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boundary conditions and some given parameters (we will specify them later). However, there 
is some difficulty in obtaining Q from the integral equation, viz. equation flMl) . So we use 
the first-order approximation, Tlj(j{r) ^ Tljfj{r + Ar) — OTmc / dr\^_^_^^ Ar. 



4. Numerical Results 

We adopt M = IOMq and Vout = lO^M in this paper. Regarding the outer bound- 
ary of the ADAF, the ion temperature should be of the same order as the virial tem- 
perature, and the electron temperature Tg should be somewhat lower than Tj because of 
the radiati on of the electrons. The angular veloc ity of the accreting fiow should be sub- 



Keplerian (INarayan. Mahadevan fc QuataertI 119981 ). So we impose the boundary conditions 



as Tj = 2 X lO^K, Te = 1 x lO^K, v/cg = 0.3. At last, by adjusting the eigenvalue of 
the problem, Nq, we can obtain the global transonic solution, i.e., a solution that can pass 
through the sonic point smoothly. 

The free parameters of our MCADAF model include a^,, cb, A, n, a, and rh, where 
m = M/MEdd with the Eddington accretion rate Msdd = 1-39 x 1O^^(M/M0)g ■ s~^ The 
first parameter describes the spin of the BH, the next three are associated with the magnetic 
field, and the last two describe the ADAF. 

Figured] shows the curve of AZd/AZh) when a* = 0.9, n = 3, A = 1 (corresponding 
to Vp = Tms), Cb = 1, a = 0.3 and m = 0.01. As the figure shows, the resistance of the 
disk is small compared with the resistance on the stretched horizon. Especially, at the outer 
boundary, the resistance of the disk is completely negligible. Since the distribution of the 
magnetic field we assumed is not smooth, there is a break at r^^. Figure [2] shows the curves 
of Q,Qf, Qh and Qp — ^ for the same parameters as Figure [H From this figure we find 
that VLf always lies between VLh and fi, which agrees with equation f|T7j) . We also find that 
the relative angular velocity of the magnetic field lines to the disk, i.e. Vlp — VL, achieves 
maximum at some radius between the inner and outer boundaries. This is a natural result 
because i^i? — ^2 is proportional to the product of two factors, AZ£)/(AZ/^ + AZ^j) and — f2 
(see equation (fTTI) ). which approach to zero at the outer and inner boundaries, respectively. 
At the innermost region, VL > VLh- This unphysical result arise because we do not use the 
exact general relativity. 

From Figure [3] we can get some ideas of the partitioning of the magnetically-extracted 
rotational energy of the BH. In this figure we show the curves of the rate of mechanical 
work due to electromagnetic torque and three kinds of heating rates per unit area. The solid 
line represents Qohm, the long dashed line for Pmw = APmw / "^T^f Ar = Hmc^, the short 



- 9 - 



dashed line for the Ohmic dissipation on the BH's stretched horizon, Qbh = AQB/f/47rrAr, 
while the dotted line for the viscous heating rate per unit area, Q^s = '^Hq^^^. As the figure 
shows, Qohm dominates over Pmw in the outer region of the disk while the latter dominates 
in the inner region. Compare Qohm + Pmw with Qbh^ we find the efficiency of extracting 
energy from the BH to the disk is very small except in the inner region. Moreover, it can be 
seen that the MC power is small in contrast to the viscous heating rate. 

The profiles of the radial Mach number, surface density S, Te and Tj, specific angular 
momentum /, and the advection factor / (= qadv/iltis + FmcI'^H)) of the accretion flow for 
different cb are shown in Figure HI The solid, dotted, and dashed lines correspond to cb = 1, 
0.5, and 0, respectively. 

From Figure H] we find that, when the MC process is present, the sonic point moves 
inward, Te decreases, S and Tj increase, while / and / decrease in the outer region and 
increase in the inner region. These effects can be understood as follows. As the spin of the 
BH is very fast {a^, = 0.9), angular momentum and energy are transferred from the BH to 
the disk. The energy flux raises the temperature of the ions. The angular momentum flux 
hinders the infalling of the accreting material. Thus the sonic point moves inward and the 
surface density increases. Since the optical depth is proportional to the surface density, the 
Compton cooling rate goes up, and consequently, the temperature of the electrons decreases. 
Compared with the case without MC process, the np component of the viscous stress tensor 
in a MCADAF around a fast rotating BH is a bit larger due to the higher pressure p in the 
outer region, so the angular momentum is transferred more efficiently and the specific angular 
momentum there is smaller. But at small radius, the increase of the angular momentum due 
to the MC process dominates over the decrease due to the viscous torque, so that the specific 
angular momentum can even increase, as can be seen from equation (pTI) . Similarly, in the 
outer region of the disk the radiative cooling rate becomes higher and / decreases, while in 
the inner region / goes up because the heating rate due to the MC process increases more 
quickly than the radiative cooling rate does. Additionally, as the magnetic field threading the 
BH becomes stronger, the influences of the MC process become more significant. However, 
as Figure [3] shows, the contribution of the MC process is of less importance, so its influences 
are small. 

Figure [5] shows the influences of the parameter A (ref. the paragraph below eq. 1). The 
solid and dashed lines are for A = 1 and 1.5, respectively. For comparison purpose the dotted 
line is shown for the case without MC process. From this figure we see that the effects of 
increasing A are similar to those of increasing cb- This is because the magnetic field in the 
region r > Vp strengthens as these two parameters increase. Although the increase of A also 
leads to the decrease of the magnetic field in the region rn < r < Vp, the MC effects are 
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small in this region, because the gravitational force there is so strong and the radial velocity 
is so high that the energy and angular momentum transferred by the MC process do not 
play any significant role. 

Figure [6] shows the effects of parameter n (ref. eq. 1). The solid and dashed lines are 
for n = 4 and 3, respectively. The lines for the case without MC process are shown with the 
dotted lines. From this figure we find that the MC effects are more significant for smaller 
n. It is because the magnetic field is weaker in the outer region (r > r^) when n is bigger, 
which can be seen from equations ([T]) and 

We calculate the influence of the MC process on the radiative efficiency of the ADAF. 
The results are shown in Figure [3 The parameters are a^, = 0.9, cb = 1.0, A = 1 and n = 3. 
The quantity M^wi denotes the critical acc retion rate of an ADAF, which is ~ a^MEdd 
(e.g. iNarayan. Mahadevan fc Quataertlll998l ). The radiative efficiency increased by the MC 
process is written as r^MC = Wcadaf — ''^adaf, where ?7mcadaf and ?7adaf are the efficiencies 
of the MCADAF and pure ADAF, respectively. From Figure [7] we find that, when a = 0.3, 
the MC process can raise the efficiency by about one percentage point, i.e., ~ 30%r7ADAF- 
But when a = 0.1 the effect of the MC process is very weak. The efficiency goes up because 
of two reasons: firstly, the MC process transports additional energy to the ADAF; secondly 
the angular momentum transported by the MC process decreases the radial velocity of the 
ADAF, and thus makes the ADAF more efficient in radiating. When a is smaller, the 
specific angular momentum of the accreting material is larger since viscosity is less efficient 
in moving angular momentum out, and consequently, the difference between the angular 
velocities of the BH and the disk is smaller. Considering equations (fT4|) . (ITS!) and (fT6|) . the 
angular momentum and energy transported by the MC process decrease. 

In all the above discussions, the spin of the BH is very large. The angular momentum 
and energy are transferred from the BH to the disk. If the BH rotates slowly, the angular 
momentum and energy may be transferred from the disk to the BH. However, since Hmc is 
proportional to a*, the MC effect is not so significant as the case when a^, = 0.9. In addition, 
there is a critical value of a^,, at which the total energy and angular momentum transmitted 
by the MC process a re zero and conseq uently rjMc = 0. In the SSD case this value is about 



a, = 0.283 forra = 3 flWang et al.l 120031 ). whereas in our MCADAF model, it is about 0.172. 
The critical value is smaller in the MCADAF case because the angular velocity of the ADAF 
is sub-Keplerian. 

The spin of the BH can even be negative, i.e., retrograde spin. From equations fll4p and 
f llSp it is easy to find that the angular momentum always flows from the ADAF to the BH 
when a^, < 0. If the resistance of the ADAF is zero, i.e. = or Qp = ^, the energy flows 
in the same direction as that of the angular momentum, as can be seen from equation f|T6|) . 
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If ^ is nonzero and big enough, Ohmic dissipation in the disk may offset the loss of energy 
that conveyed to the BH. According to equation f|T6|) . the critical condition is APmc = 0, 
or equivalently, ^ + Phd = 0. Our calculations show that when a^, < 0, ^ is so small that 
^ + Phd < holds for almost all cases, and the net energy flows from the ADAF to the BH. 
If the BH rotates rapidly, e.g. a^, = —0.9, the effects of the MC process will be negative 
compared to the case of a* = 0.9: the radial velocity and the temperature of the electrons 
increase, the sonic point moves outward, the temperatures of the ions decrease, etc. 



5. Summary and Discussion 

In this paper we investigate the influence of the magnetic coupling process on the dy- 
namics of the ADAF. The effects of the MC process on the basic equations of the accretion 
flow is simplified as a source of angular momentum and energy. The angular momentum 
and energy fluxes are derived with the equivalent circuit approach. We find that when the 
BH rotates fast (e.g., a^, = 0.9) and when the viscous parameter a is large, a = 0.3, for 
reasonable magnetic field, the MC process can mildly affect the dynamics of the ADAF, 
increasing the ion temperature and the density of the accretion flow and decreasing the elec- 
tron temperature and the radial velocity. The MC process can also raise the effciency of the 
ADAF by ~ 30%. But if a = 0.1 or smaller, the influences of the MC process are much 
weaker and can be neglected. 

In the above calculat i ons th e strength of the magnetic field is estimated following 



Moderski. Sikora fc Lasotal (jl997l ). which is equivalent to assume < 1- Obviously un- 
certainties exist in the above estimation. On the one hand, recent MHD simulations show 
that the magnetic field strength near the horizon can be very high, almost four times as large 



as the the equipartition value (iMcKinneyl l2005l ) . which corresponding roughly to cb ~ 2. 



On the other hand, our calculation requires that there exists an upper limit to the value of 
Cb- This is because if cb were too large, the transferred angular momentum from the BH 
to the accretion flow would be so significant that the accretion could not proceed due to the 
strong centrifugal force. We find that the highest value of cb depends on the accretion rate 
for given a^,. It can be ~ 5 if the accretion rate is very low but ~ 1 if the accretion rate is as 
high as m > 0.1. Considering the above two limitations on cb, the increased efficiency due 
to the MC process ?7mc(= ''^mcadaf — ''7adaf) can be as high as 10%. 

Observations of the hard state of BH X-ray binaries sometimes show luminosities as high 
as ~ 10 — 30%LEdd, which are much higher than A VoLEdd, the highe st luminosity that an 



ADAF surrounding a Schwarzschild BH can produce (lEsin et al.lll997l ). When the accretion 



rate is higher than the critical rate of an ADAF, the accretion flow enters into the regime 
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of the Luminous Hot Accretion Flow (LHAF) flYuanlboOlh . lYuan et all (120071 ) found that 
the highest luminosity an LHAF surrounding a non-rotating BH (with M = O.SMEdd) can 
produce was ~ 8%LEdd5 which is still too low. They speculated that when the spin of the BH 
and the MC process were taken into account, the highest luminosity an LHAF could produce 
would possibly be high enough to explain the observed high Lx- Now, taking into account 
the black hole spin and the MC process, we recalculate the maximum luminosity an LHAF 
(with rfi = 0.3) can produce, with reasonable parameters such as a^: = 0.9, A = 1, ~ 1 
(note this is the largest possible value when m = 0.3) and n = 3. We find the highest 
luminosity is ~ 14%LEdd- The increase due to the MC process is ~ 1.0%LEdd while that 
due to the BH spin is ~ 5.1%LEdd- So the MC process seems not so helpful to increase the 
highest luminosity an LHAF can produce to explain the observed highest L^. 

In this paper we assume that the closed field extends to the outer boundary of the 
ADAF according to a power law. As a matter of fact the magnetic connection b e tween 
the BH and the disk can be maintained only within a limited radius. I Wang et al.l (120041 ) 
discussed the constraint of the screw instability to the MC region of a SSD based on the 
Kruskal-Shafranov criteri on: the screw ins t ability will occ ur, if the magnetic field line turns 
around itself about once (jKadomtsevlll966l : lBatemarull978l ). It turns out that the MC region 
is limited within some critical radii on the SSD. By numerically solving the Grad-Shafranov 
equation, t he ma in differential equation that describes the structure of the magnetosphere, 
Uzdenskyi ( 120051 ) argued that for a rapidly-rotating BH the field lines are frame-dragged 
by the BH so much, and the toroidal magnetic field becomes so strong that the magnetic 
connection between the BH and the disk cannot be maintained over a large range of radii 
on the disk. We shall address this issue in the context of ADAF in our further work. 
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Fig. 1. — Profile of ^(= /S.Zd/ ^Zh) on the disk. The parameters are a* = 0.9, n — Z^X — l 
(corresponding to Vp — Vms), cb — 1, Oi — 0.3 and m — 0.01. 
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Fig. 2. — Curves of the angular velocities. The solid, dashed and dotted lines correspond 
to Qp, Q and Qh, respectively. The thick solid line shows flp — fl. The parameters are the 
same as Figure 1. 
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Fig. 3. — Curves of the rate of mechanical work due to electromagnetic torque and the 
heating rates per unit area. The solid line represents the rate of Ohmic heating in the disk 
Qohm, the long dashed line shows the power of the electromagnetic torque on the disk Pmw, 
the short dashed line shows the Ohmic dissipation on the BH's stretched horizon Qbh- As 
comparison the curve of viscous heating rate per unit area is shown in dotted line. The 
parameters are the same as Figure 1. 
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Fig. 4. — The profiles of the radial Mach number, specific angular momentum /, electron and 
ion temperatures Tg and Tj, surface density S (in unit of Sq = M/M), and the advection 
factor / for different cb- Other parameters are the same as Figure 1. The solid, dotted, and 
dashed lines correspond to cb = 1, 0.5, and 0, respectively. The thick solid line in the top 
right panel shows the Keplerian angular momentum. 
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Fig. 5. — Same with Figure 1, but for different A when a* = 0.9, cb = 1, a = 0.3, and 
771 = 0.01. The sohd and dashed hues correspond to A = 1 and 1.5, respectively. The dotted 
hne corresponds to the case without MC process. 
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Fig. 6. — Same with Figure 1, but for different n when a* = 0.9, cb = 1, a = 0.3, A = f .5, 
and 771 = 0.01. The sohd and dashed hues correspond to n = 4 and 3, respectively. The 
dotted line corresponds to the case without MC process. 
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Fig. 7. — The radiative efficiency of the accretion flow as a function of the accretion rate. 
Here Merit = a^MEdd- The thin sohd fines correspond to the results of the present MCADAF 
model (a* — 0.9, = 1.0, A = 1.0 and n = 3), the dashed lines indicate the results without 

MC process, the thick solid lines show the efficiency improved by the MC process {f]Mc), 
and the dotted lines show the results of a Schwarzschild BH without considering the MC 
process. In the upper panel a = 0.3 and in the lower panel a = 0.1. 



